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ABSTRACT 


The Barnes PRT-5 Infrared Radiometer and the L-Band (1.42 GHz) 
Channel of the Multifrequency Microwave Radiometer (MFMR) are used to 
survey the distribution of surface water temperature and salinity. 

These remote sensors were flown repetitively in November 1971 over the 
outflow of the Mississippi River into the Gulf of Mexico. Data reduc- 
tion parameters were determined through the use of flight data obtained 
over a known water area. With these parameters, the measured infrared 
and microwave radiances were analyzed in terms of the surface tempera- 
ture and salinity. 
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INTRODUCTION 


Surface water salinity and temperature measurements are useful to 
studies of coastal processes such as hay circulation and river outflow. 
Marine organisms are affected by the salinity and temperature of their 
water habitat. Marine ecologists and biologists thus need the capa- 
bility for conducting surveys of salinity and temperature. Also, fish 
forecasting by the National Fisheries Service of the National Oceanic 
and Atmospheric Administration (NOAA) requires information on the tem- 
perature and salinity of breeding grounds of marine organisms . 

Past investigations have shown that the microwave emission of 
water is affected significantly by changes in the salinity of the water 
surface layer for microwave frequencies less than 3 GHz(Rfs. 1-3). 

Thus, the possibility existed for making surveys of the water surface 
salinity from aircraft through the use of an L-band microwave radi- 
ometer. A Multifrequency Microwave Radiometer (MFMR) built by Aerojet 
General Corporation was installed on the NASA 927 (Lockheed P3A aircraft) 
in 1969 . One channel of the MFMR operates in the L-band at 1.42 GHz. 

The purpose of this paper is to report on the development of an analysis 
technique whereby the distribution of surface water temperature and 
salinity along the flight line may be calculated by computer from MFMR 
and infrared radiometry (PRT-5) data taken along the flight track of 
the NASA 927 aircraft . 

In order to recover the salinity of the water surface from MFMR 
measurements , one must consider carefully the effects of a number of 
other parameters, such as, variations in water surface temperature, 
absorption and emission by the intervening atmosphere, absorption and 
emission by the MFMR radome antenna, and waveguide, and departures in 
the calibration constants from the values determined before the air- 
craft flight. These considerations are discussed in Section II below. 

A computer program was developed to analyze the MFMR and PRT-5 
data and to produce profiles of water surface salinity and temperature 
along the aircraft flight track. The development of the computer pro- 
gram is discussed in Section III below. 

The computer program was used with flight data acquired during a 
flight of the NASA 927 over the outflow of the Mississippi River into 
the Gulf of Mexico in November 1971. The results of the application 
of the program are given in Section IV below. 


Conclusions and recommendations are presented in Section V below. 
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THEORETICAL DEVELOPMENTS 


EMISSION AND REFLECTION OF WATER 


In the determination of salinity from measurements of the microwave 
emission of water, one needs a computationally simple, mathematical 
model that may be used to predict the microwave emissive and reflective 
properties of water, given the observable conditions of the water and 
trial values of salinity. It is known that the microwave emission of 
water is affected, in general, by changes in water temperature, water 
salinity, water surface roughness (wind waves and swell) , and by the 
presence of other materials on the surface, such as sea foam and oil. 
(Refs. 1-7) 

Mathematical models exist for predicting the microwave properties 
of sea foam and oil on water. These models require, however, a know- 
ledge of the thicknesses of foam and oil, of the dielectric constant of 
the oil, and of the air-to-water ratio of the foam. These parameters 
are not available through the use of remotely sensed data. On the other 
hand, visual or photographic data may be used to determine the presence 
or non-presence of oil or foam on the water surface. Thus, in the 
mathematical model, it is assumed that there is no foam or oil on the 
water surface. 

Mathematical models exist also for predicting the effect of water 
surface roughness on the microwave emission of water. These models 
require a knowledge of the spectrum of water surface slopes which, in 
turn, may be inferred from a knowledge of the fetch, duration, and speed 
of the wind over the water. Again, these parameters are not available 
from remotely sensed data. Studies (Refs. 4 and 6) have shown, however, 
that the effect of the water surface roughness is small for viewing 
angles near nadir and for vertical polarization. In the mathematical 
model developed for use in this paper, it is assumed that the effects 
of water surface roughness can be ignored. 

With the elimination of the above effects, only the effects of 
temperature and salinity remain. In order to model the microwave prop- 
erties of water, it is assumed that the water is a homogeneous, semi- 
infinite dielectric bounded by a plane surface and overlaid by air. 

Since the water is an absorbing and emitting medium of large depth, the 
radiation impinging upon the air-water interface from below is the same 
as that emitted by a black-body at the temperature of the water. That 
is, the brightness temperature of the radiation just under the air-water 
interface is equal to the absolute themometric temperature of the water. 
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At the air-water interface, a fraction, r , of the radiation 
impinging from below is reflected back into the water mass. The rest 
of the radiation is transferred through the air-sea interface and 
emerges as emitted radiation. Thus, the brightness temperature of the 
emitted radiation, T^ (°K), is 

T = (1 - r) T (1) 

e w 


where T^ is the water temperature (°K). 

The reflectivity, r , for a flat semi-infinite dielectric of 
dielectric constant, e'-je", is given as (Ref. 2) 
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where \p is the zenith angle of observation and the subscripts v and 
h refer to vertical and horizontal polarization, respectively. 

The microwave radiation directed downward impinging upon the air- 
water interface is characterized by a sky brightness temperature, T sk y 
(°K), which is the sum of the transmitted cosmic background bright- 
ness and the katabatic emission by the atmosphere. 

A portion of the sky radiation is reflected from the air-water 
interface toward the observer, such that the total radiation leaving 
the air-water interface and traveling toward the observer has a bright- 
ness temperature of 
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In order to compute the reflectivity for a given angle of observa- 
tion and polarization, one must be able to specify the dielectric 
constant, (e'-je") of the water. 

Studies (Refs. 8 and 9) have shown that the values of e' and e" 
for water with dissolved salt (NaCiL) may be predicted for the microwave 
spectrum through the use of Debye model. That is. 
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where £ is 4.9, £ , T, and a. are functions of water temperature 

0° S 1 r 

and salinity only, and 0) is the circular frequency , [2 it V, where V 
is the electromagnetic frequency (Hz)]. 

The values of £ s , T, and Op for specific water temperatures and 
salinities as given in Refs. 8 and 9 are given in Table 1. In order 

to predict the values of e' and e" for any water temperature and 

salinity, it is necessary to interpolate between the values given in 
Table 1. Interpolation formulas exist for this purpose and are indi- 
cated in Table 2 (Ref. 10 ). 

With these interpolation formulas and the use of Equations 2-6, 
one has a mathematical model that may be used to predict the emissive 
and reflective properties of water at 1.42 GHz for any given water tem- 
perature, water salinity, angle of viewing and polarization. This model 

was used to compute the microwave emission for vertical and horizontal 
polarization, for viewing angles of 0, 30, and 55 degrees, for water 
salinities of 0, 5, 10, 15, 20, 25, 30, and 35 0/00, and for water 
temperatures ranging from 0 to 30°C. The results of the calculations 
are shown in Figures 1-5- 

These figures show that salinity has the maximum effect on micro- 
wave emission for angles of viewing near 55 deg, for vertical polariza- 
tion, for warm water, and for salinities in the range from 10 to 
35 0/00. 
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EMISSION AND ABSORPTION BY THE ATMOSPHERE 


A mathematical model that may he used to predict the effect of the 
atmosphere on the transfer of microwave radiation is needed to allow 
the proper consideration of this effect. For a given aircraft flight, 
soundings of atmospheric temperature and humidity are available usually 
through the use of radiosondes released at National Weather Service sta- 
tions. At 1.42 GHz, the atmosphere is practically transparent for all 
conditions including clouds and rain (Ref. 9)- In the case of rain, 
however, the surface water would be disturbed and diluted by the rain 
water. Also the aircraft radome would be covered by liquid water and 
would create large and unpredictable losses. For these reasons, rainy 
weather conditions must be ruled out for the measurement of salinity by 
the MFMR. The emission and transmission of the atmosphere may be com- 
puted from the knowledge of the vertical distribution of temperature 
and humidity in the atmosphere (Ref. 10 ). 

In addition, the effects of refraction and a spherical geometry may 
be easily included in the calculation. (Ref. ll). As a result, the 
relationship between the brightness temperature of the radiation 
impinging upon the MFMR radome, T^ , and T Q , the brightness temper- 
ature of the radiation leaving the air-water interface, is 


T, = T t + T (T) 
b o atm 

where t is the transmission coefficient for the intervening atmosphere 
and T atm is the brightness temperature of the radiation emitted by the 
intervening atmosphere toward the observer. 


TRANSMISSION OF SENSOR COMPONENTS 


If one neglects the radiation received in the side lobes of the 
L-band radiometer of the MFMR, then the equivalent radiometric tempera- 
ture of the radiant power received by the radiometer is given by 
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where 


and 


L is the loss factor for the radome, 

r 

L is the loss factor for the antenna, 
a 

L is the loss factor for the waveguide, 
w 

T is the thermometric temperature of the radome, 
rd 

T is the thermometric temperature of the antenna 

3 . 

T is the thermometric temperature of the waveguide, 
w 


5 


CALIBRATION EQUATION 


The microwave radiometer is a linear device such that the output, 
C ( counts J* , is given by 


C 



( 9 ) 


where C is the output of the radiometer when the power emitted by 
the reference load is being monitored, (counts), 

T is the thermometric temperature and equivalent radiometric 
k temperature of the reference load, (°K) , 

and S is the response of the instrument or scale factor (°K per count). 

In the MFMR , the scale factor may be determined before or after a 
data run by injecting a known amount of argon noise that has an equiva- 
lent noise temperature, T n (°K), and observing the change in output, 

AC (counts). In this case. 


S 



SUMMARY OF MATHEMATICAL MODELS 


( 10 ) 


Mathematical models were developed above to predict the effects of 
all measurable parameters and salinity on the output of the MFMR. These 
models are valid only for the following conditions . 

*The output of the MFMR is voltage, however, it is recorded as 
PCM counts . 
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1. Smooth water surfaces free from foam and oil 

2. Non-raining atmospheres 

3. Matched radiometer systems 

4. Linear radiometer syst ems 


PROGRAM DEVELOPMENT 


DATA AVAILABLE FOR ANALYSIS 


The data available for analysis consists of the following: 

1. One-second averages of the thermometric temperature of the radome, 
the antenna, and the waveguide of the MFMR and of the equivalent 
radiometric temperature of the radiation measured by the L-band 
radiometer as determined from the calibration equations (Eqs. 9 
and 10). 

2. Averages taken over 0.792 sec. of the equivalent radiometric temper- 
ature of the radiation measured by the PRT-5 infrared radiometer. 

3. Radiosonde measurements of atmospheric temperature and humidity from 
the surface to the 10 mb level (approximately 31 km). 

b. Flight logs for the aircraft flight. 

5. Location (latitude, longitude, and altitude), attitude (heading, 
drift, roll, pitch), and speed of aircraft sampled once each ten 
seconds . 

6. Surface measurements of salinity and temperature at selected points 
under the aircraft surface track. 

7. Photographic coverage (10% overlap) along flight line. 

8. Scale factors, baseline counts, and baseline temperature used to 
reduce raw MFMR data to produce data in 1 above. 

9. Equivalent argon noise temperature, radome loss factor, and wave- 
guide loss factor. 
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ADEQUACY OF THE AVAILABLE DATA AND MATHEMATICAL MODELS 


In the above it has been noted that some minor effects such as 
surface roughness and sidelobe contributions have been ignored in the 
development of the mathematical models. In addition, there is some 
evidence that there are systematic errors in the values of injected 
argon noise and of the antenna loss factor for the MFMR. All of these 
factors could create a significant bias in the data if not properly 
treated. In order to insure consistency in the analysis and to elimi- 
nate the effects of systematic errors, a new set of calibration con- 
stants (scale factor and antenna loss factor) are obtained through the 
use of flight data and the mathematical models. Two portions of the 
flight data are taken — the first over a stretch of water of known 
temperature and salinity and the second when the radiometrically cold 
sky is being viewed. 


DETERMINATION OF THE NEW CALIBRATION CONSTANTS FOR THE MFMR 


If equations (8) and (9) are combined and rearranged, one obtains 
the expression. 


T = L L L [T . - S (C . - C)1 
b raw 1 b£ v b£ /J 
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The actual output of the radiometer , C , may be obtained from the 
given values of radiometric temperature through the use of equiations 
(9) and the old calibration constants. 

If T^ is the brightness temperature corresponding to the radi- 
ometer reading 5 ^w > obtained when the known water target is being 
viewed, and T^g is the brightness temperature of the sky corresponding 
to the radiometer reading, C s , then from equation (ll). 
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Inserting equation (13) into (ll) and solving for L , one 
obtains the expression 
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In equations (l4)~ (l6), one may use either C w and Tj-, w or C s 
and Tjjg to determine L a . With the knowledge of L a , one may 
obtain the value of the scale factor, S , from equation (13). 

The effect of the above determinations is to force the MFMR measure- 
ments to agree with the mathematical model at two points, namely, for a 
selected water point and for a cold sky point . Since all of the radio- 
metric measurements are made for water conditions close to the selected 
water point, there should be consistency in the values of the derived 
salinities . 


COMPUTATIONAL PROCEDURE FOR OBTAINING SALINITY AND 
TEMPERATURE FROM MFMR AND PRT-5 DATA 


After receipt of all of the flight data and surface and atmospheric 
measurement data obtained in support of the aircraft flight, the fol- 
lowing procedure was followed: 

1. A homogeneous area of water with known water temperature and 
salinity was selected. 
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2. From an inspection of the surface measurement and PRT-5 tabulations, 
the difference between the two temperatures was determined for the 
selected water area. In general, there will be a fixed difference 
due to the cool skin of the water, the nonunity value of the emis- 
sivity of water, and the absorption and emission of the intervening 
atmosphere . 

3. The radiation properties of the atmosphere were calculated through 
the use of the Microwave Atmospheric and Oceanic Data Analysis 
Program (MAODAP). The MAODAP uses radiosonde data taken during the 
flight (Ref. 12). 

4. The salinity and temperature analysis program (STANL) was then 
employed. In the STANL, the following operations are performed: 

a. Controls cards are read. 

b. Data are^read from PRT-5 and MFMR data tapes supplied by the 
Computing and Analysis Division (CAAD) of NASA Manned Space- 
craft Center. 

c. Eleven-second running averages are computed for needed MFMR 
and PRT-5 parameters . 

d. The new calibration constants are determined by use of the 
flight data and the Eqs . 13-l6. 

e. The salinity and temperature are determined from the MFMR 
and PRT-5 data through the following procedure: 

(1) First, the microwave radiometric temperatures are 
corrected for the change in scale factor, that is 

(T - T ) 

T' = T + S' — ~ — ~ — (IT) 

r bJo S 

(2) Second, the PRT-5 infrared radiometric temperature are 
used with the correction factor in 2 above to determine 
the water temperature (°K). 

(3) Last, the value of salinity is determined through itera- 
tion such that the predicted radiometric temperatures 
agrees with the corrected radiometric temperatures 
determined in (l) above. 
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PROGRAM APPLICATION 


MISSION 190, FLIGHT 3 


Mission 190, Flight 3 was flown over the outflow of the Mississippi 
River into the Gulf of Mexico (Site 128) on November 11, 1971. Nine 
data runs were made over Line 1 (see Figure 6 ) at approximately 242 m 
(860 ft.) above the water surface. The first eight data runs were made 
with the MFMR antennas pointing 4.8 deg from nadir in vertical polari- 
zation. The last run was made with the antennas pointing 45 deg from 
the zenith in vertical polarization. The first data run began at the 
ENE end of Line 1 at 0935 CST; the last data run ended at the ENE end 
of Line 1 at 1803 CST. Two KA-62 cameras, the Barnes PRT-5 Infrared 
Radiometer, the MFMR Boresite Camera, the MFMR, and the navigation 
equipment were used on the NASA 927 (Lockheed P3A) aircraft during the 
flight, and all sensors functioned normally. 

Flight support consisted of a radiosonde released at Boothville, 

La. (approximately 25 n. miles north of the aircraft flight line), and 
a surface survey along the aircraft flight line. The surface survey 
was made through the use of a 65 foot crewboat leased at Venice, La. 

The surface measurements consisted of in situ measurements of salinity 
and temperature with a Beckman In-Situ Salinometer. The measurements 
were taken at stations approximately every nautical mile along the 
flight line. In addition, water samples were taken approximately every 
ten stations. As originally planned, the surface survey was to take 
place during a time period centered upon the time interval of the air- 
craft flight. Unfortunately , the crewboat ran into the dock the evening 
before the day of the aircraft flight and had to undergo repairs before 
leaving during the morning of the flight. The accident did not come to 
the attention of the mission manager since the aircraft crew and sensor 
operators had left for the airport before the surface survey team 
arrived at the Gulf Docks in Venice to leave on the crewboat. As a 
result, the surface survey boat arrived at S. Pass and began the survey 
at 1157 CST just before the aircraft survey was being completed. The 
surface survey proceeded toward the WSW end of Line 1 and returned to 
S. Pass at l64l CST. 

In general, the atmosphere was clear and dry, and surface winds were 
from the ENE at 5 knots. The radiosonde was released at 1715 CST. The 
results of the radiosonde sounding are shown in Table 3. 
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DATA PREPARATION 


The water area at the Southwest Pass outflow was selected as the 
calibration point for the water data^ The surface measurements of water 
surface temperature at this point were compared with the PRT-5 radio- 
metric temperatures , and the radiometric temperature was found to be 
3.5 °C lower than the bucket temperature of the water surface. 

The MAODAP was used to calculate the microwave parameters of interest 
for the subsequent analysis. The parameters of interest from the MAODAP 
program are listed in Table 4. 

The analysis program (STANL) developed for the study was then used 
to determine the new calibration constants and to produce a plot of sur- 
face water salinity and temperature along the flight track. 


RESULTS 


The results of the determination of new calibration constants for 
the MFMR are listed in Table 5 T A significant difference was noted 
between the old and new values of antenna loss factor, scale factor and 
injected noise. These differences appear to be the result of a bias in 
the measurement of injected noise as made by the hot/cold load labora- 
tory calibration procedure. 

The plot of surface water salinity and temperature along the flight 
track during run 8 is shown in Figure 7. The symbols in the Figure 
indicate the values of surface temperature and salinity obtained by 
the surface survey that was conducted shortly after the aircraft flight. 
In general, there is good agreement. The times of passage over the 
river mouths are indicated by the symbol M and the time of passage 
over surface foam lines as seen in the photography are indicated by 
the symbol F . 

These data confirm the usefulness of L-band measurements of water 
microwave emission in determining the surface salinity of water. 
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TABLE I. -THE DEBYE PARAMETERS OF AQUEOUS SODIUM CHLORIDE 
FOR SPECIFIC VALUES OF TEMPERATURE AND SALINITY 
(AFTER SAXTON AND LANE, 1952 5 SAXTON, 1952) 
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Debye Parameter 
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0.0 
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TABLE II. -REGRESSION COEFFICIENTS FOR THE DEBYE PARAMETERS 
OF AQUEOUS SODIUM CHLORIDE 


Independent 

variable 


Regression coefficients 


e 

T x 10 12 (sec) 

0 . (U m 1 ) 


s 

t 
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constant 

0.88195 E< 02 

T 

sea 

-0.40349 E-00 

S 

-0.43917 E-00 

T S 

sea 

0.43269 E-02 

T 2 

sea 

0.65924 E-03 

s 2 

0.16738 E-02 

T 2 S 
sea 

- 0.92286 E-05 

2 

T S 

sea 

-0.42856 E-04 

T 2 S 2 
sea 

eX P (T sea } 

0.44410 E-07 


0.19390 E 02 


-0.68020 E-00 


-0.11370 E-00 

0.87483 E-01 

0.58629 E-02 

0.45802 E-02 

0.95865 E-02 


0.11417 E-02 

00.25662 E-04 

-0.87596 E-04 

-0.16914 E-04 

-0.54577 E-04 

-0.37158 E-04 

0.82521 E-06 

0.39288 E-06 

-0.65303 E-17 
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TABLE III. -RADIOSONDE DATA OBTAINED AT BOOTHVILLE , LA., AT 
1715 CST , 11 NOVEMBER 1971, IN SUPPORT OF 
MISSION 190, FLIGHT 3, SITE 128 


Pressure 

(mb) 

Temperature 

(°c) 

Dew-Point Depression 

(°c) 

1023.5 

17.2 

5.3 

1000. 

17.3 

5.6 

953. 

13.5 

5.1 

939. 

13.8 

27.0 

850. 

12.2 

28.6 

776. 

9.1+ 

28.0 

500. 

-13.1* 

21.9 

1 * 00 . 

-26.9 

18.3 

3lU. 

-39.4 

15.8 

200. 

-61.7 

Dry 

171. 

-65.1* 

Dry 

159. 

- 63.8 

Dry 

128 . 

-68.1* 

Dry 

121. 

-66.1 

Dry 

100. 

-70.6 

Dry 

95. 

-73.0 

Dry 

50. 

- 62.5 

Dry 

1*1. 

-55.3 

Dry 

23. 

-53.5 

Dry 

17. 

-1*6.2 

Dry 

13. 

-1*7.6 

Dry 

10. 

-1*0.5 

Dry 



TABLE IV. -LIST OF PARAMETERS FROM AMODAP NEEDED 
FOR THE ANALYSIS PROGRAM 




Parameter 

Value 


Temperature of selected water area 

22.2 °C 

Salinity of selected water area 

8.0 0/00 

Equivalent radiometric temperature of sky 
(at 1+5 deg zenith angle) 

7.06 °K 

Atmospheric transmission factor for intervening 
atmosphere 

0.99982 

Emissive brightness temperature of intervening 
atmosphere 

0.1 °K 

Sky brightness temperature at surface for the 
viewing angle of the MFMR 

5.15 °K 

Zenith angle of viewing 

4.8 deg 
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TABLE V . -CALIBRATION CONSTANTS FOR THE MFMR 


Parameter 

Old 

Value 

New 

Value 

Radome loss factor 

1.10786 

1.10786 

Antenna loss factor 

1.600 

1.509 

Waveguide loss factor 

1.0423 

1.0423 

Scale factor 

0.5875 °K/cnt 

0.49588 °K/cnt 

Injected noise 

133.0 °K 

112.3 °K 



Emissive Brightness Temperature of a Calm 
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Temperature of Sea Water, T sea (C) 

Pig. i. Microwave emission of a calm sea: v — 1.42 GHz, 

0 =0 deg, vertical or horizontal polarization. 






o 5 10 15 20 25 30 

Temperature of Sea Water, T sea (C) 

Fig. 3. Microwave emission of a calm sea: v = 1.42 GHz, 

9 = 55 deg, vertical polarization. 
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FLIGHT LINE FLOWN DURING MISSION 190, 
FLIGHT 3, NOVEMBER 11, 1971 



MISSISSIPPI DELTA SITE 128 


Fig 6. Flight line flown during mission 190, flight 3, 
site 128, November 11, 1971. 
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Fig 7. Results of mission 190, flight 3, run 8 
over site 128, line 1. 




